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1. Introduction 

The structure of hadrons, such as size, shape and charge distribution can be probed by the 
electromagnetic form factors. The CI baryon, consisting of three valence strange quarks, is sig- 
nificantly more stable than other members of the baryon decuplet, such as the A, with a life-time 
on the order of 10~ 10 s. This fact makes the calculation of its electromagnetic form factors par- 
ticularly interesting since they are accessible to experimental measurements. Its magnetic dipole 
moment is measured to very good accuracy unlike those of the other decuplet baryons. A value of 
/ift- = —2.02(5) is given in the PDG [p in units of nuclear magnetons (/%). Within lattice QCD 
one can directly compute hadron form factors starting from the fundamental theory of the strong 
interactions. Furthermore, higher order multipole moments, not detectable by current experimental 
setups, are accessible to lattice methods. Higher-order moments such as the electric quadrupole 
are essential in the determination of the deformation of a hadron state. In this work, we calculate 
the electromagnetic form factors of the Q. baryon using, for the first time, dynamical domain-wall 
fermion configurations. For the calculation we use the fixed-sink approach which enables us to 
calculate the form factors for all values and directions of the momentum transfer q concurrently. 
The main advantages of this approach is that it allows an increased statistical precision, while at 
the same time it provides the full Q 2 dependence. Moreover, by constructing optimized sources 
we calculate the two dominant form factors accurately. An appropriately defined source for the 
subdominant electric quadrupole form factor is constructed and tested [0]. The aim is to obtain an 
accurate determination of the quadrupole moment, at a price, of course, of additional sequential 
inversions. 



2. Electromagnetic form factors of the Q. baryon 



The on-shell £1 matrix element of the electromagnetic current j^ M , can be decomposed in 



terms of four independent Lorentz covariant vertex functions, a\{q ), a 2 {q ), c\{ 



which depend only on the squared momentum transfer q 2 
spacetime these are given by [||] 
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The rest mass and the energy of the particle are denoted by niQ and The initial (final) four- 
momentum and spin-projection are given by Pi(pf) and Si(sf). In addition, every vector com- 
ponent of the spin- 1 Rarita-Schwinger vector-spinor u G satisfies the Dirac equation, {p^J^ — 
m£i)u° \p,s) = 0, along with the auxiliary conditions: y a u c {p : s) = and p a u a (p,s) = 0. Further- 
more, the covariant vertex functions are linearly related to the experimentally measured electric 
GEo{q 2 ), GEi{q 2 ) and magnetic Gmi {q 2 ) and GM?>{q 2 ) multipole form factors [[|, |[|. 
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3. Lattice techniques 

We employ gauge configurations generated by the RBC-UKQCD joint collaborations using 
Nf = 2 + 1 dynamical domain-wall fermions (DWF) and the Iwasaki gauge-action corresponding 
to a pion mass of about 330 MeV. The simulation is carried out on a lattice of size 24 3 x 64 with 
a lattice spacing of 0. 1 14(2) fm. For this pion mass and finite lattice volume the Q. is stable. 
The lattice spacing a, the light u and d and the strange quark mass were fixed by an iterative 
procedure using the £2 , the pion and the kaon masses 0]. For the present calculations we use 200 
well separated dynamical domain-wall fermion gauge configurations It is known that the chiral 
symmetry breaking falls exponentially with the length of the fifth dimension of the DWF-action. 
The value Ns = 16 used here is adequate to keep the residual mass sufficiently small. 

3.1 Interpolating fields 

In order to calculate the Q. matrix element we need to evaluate the appropriate two- and 
three-point correlation functions. An interpolating field operator with the quantum numbers of the 
£2~ hyperon is given by 



(3.1) 



where C = 74/2 is the charge-conjugation matrix and o is the vector index of the spin-3/2 spinor. 
In order to ensure ground state dominance for 
the shortest time evolution we perform a gauge 
invariant Gaussian smearing, as described in 
Refs. [Bl H], on the quark fields s that enter in the 
interpolating field: 



qp(t,x)=Y^ + aH{x,y-U)] n qp{t,y), 



H(x,y;U)=Y t (L^(x,0%-£ + ^(*- A,0%+#) 

where q (q) is the smeared (local) fermion field. 
The links U^(x, t) entering the hopping matrix H 
are APE-smeared gauge fields. 
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Figure 1 : The Q. effective mass and the fit to 
the plateau plotted against time separation. 



For the lattice spacing considered here we have used the Gaussian smearing parameters a = 
5.026 and n = 40, which have been optimized for the nucleon state. In Fig. 1 we show the CI 
effective mass, calculated from the two-point function ratio defined by ra^y(?) = — log[G(? + 
1,0)/G(?,0)]. It displays a nice plateau yielding mp, = 1.76(2) GeV. This value is 5% higher 
than the experimental one. This may reflect a slightly larger value for the strange quark mass as 
compared to the physical one. 



3.2 Two- and three-point Correlation functions 

The electromagnetic form factors are extracted by constructing appropriate combinations of 
two- and three-point correlation functions. The corresponding lattice correlation functions are 
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given by 

3 



G(t f ,q) =I£e- S Hr 4 ^(^(x / )/ ya (0)), (3.2) 
Gi\(^^)=£y*^^ (3.3) 

XfX 

where F 4 = + y 4 ) and = jl^y 5 ]^, with the Dirac y-matrices taken in Euclidean spacetime. 
We calculate the above three-point correlation function in a frame where the final Q. state is at 
rest, pf = 0, while the operator j EM is inserted at time t carrying a momentum q = —pi. 

The leading time dependence and unknown overlaps of the Q. state with the initial state Jq\0) 
in the three-point correlation function can be canceled out by forming appropriate ratios of the 
three-point function and two-point functions. A particularly suitable ratio is defined by 



R n (r - 1) = G a ^(r v ,q,t) G kk (r\p h t f -t)G kk (r\o,t)G kk (r\o,t f ) 

at[ ,<? ' G kk (T*At f ) \ G kk {T\Q,t f -t)G kk {T\piS)G a {T\p h t f y 

where a summation over the repeated indices k (k = 1,2,3) is assumed. For large Euclidean time 
separations this ratio becomes time-independent (plateau region) 

R a \(T,q,t) 'i-'^ 1 n a \(F,q) = V Tr [TK^ipf^'h^ipi)] , (3.5) 
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It is understood that the trace acts on spinor-space, while the Euclideanized version of the Rarita- 
Schwinger spin sum is given by 
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The electromagnetic form factors are extracted by fitting to the plateau region of a set of 
specially chosen combinations of three-point functions given below 



£ nf k (r\q) = K V G EQ (Q 2 ) +K 2 G E2 (Q 2 
k=\ 

£ Sjk,nJ l k (r 4 ,q)=K 3 G Ml (Q 2 ), 

j,k,i=i 

^e jkl U J \^,q)=K 4 G E2 (Q 2 ). 
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(3.8) 
(3-9) 
(3.10) 



The continuum kinematical coefficients Ki (i = 1,2,3,4), have been calculated in Ref. ^ and are 
functions of the £2 mass, the energy Eq, the space-time index p. and the momentum q. Further- 
more, it is readily seen that the subdominant electric quadrupole form factor G E2 is isolated by 
the combination provided by Eq. ( 3.10p . In addition, these expressions are such that all possible 
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directions of jj, and momentum q contribute in a symmetric fashion at a given momentum transfer 
Q 2 . In other words, the optimal combinations of IT^ t (T,q) employed are those that maximize the 
number of nonzero q contributions in a lattice rotationally invariant way [§J. 

We calculate only the connected contributions to the three-point function by performing se- 
quential inversions through the sink. This means that we need to fix the quantum numbers of the 
initial and final states and explains why we consider optimal combinations. Otherwise, one would 
have to prepare a new set of sequential inversions for every choice of vector and Dirac indices, 
which is prohibited in a lattice computation given the fact that we have to consider an overall of 
CJxj3xTxa = 256, combinations by looking the index structure of Eq. (3.3). Having the matrix 
element for all the different directions of q and for all four directions /I of the current, we form 
an over-constrained system of linear equations (in terms of form factors), which is solved by em- 
ploying a singular value decomposition analysis. This procedure yields the electromagnetic form 
factors Geo, Gmi an d Gei- The statistical errors are found by a jack-knife analysis, which takes 
care of the correlations of lattice measurements of the ratios. 

As already mentioned the three-point function of the connected part is calculated by perform- 
ing sequential inversions through the sink. This requires to fix the temporal source-sink separation. 
In order to determine the smallest time separation that is still sufficiently large to damp the excited 
state contributions we perform the calculation at two values of the sink-source separation, namely 
tf/a = 8 and tf/a = 10. We compare in Fig. [2| the results for the plateaus Uff X (T,q), for a few 
selected directions of the current and for low momentum q values for these two sink-source time 
separations. As can be seen, the plateaus values at tf/a = 10 are consistent with the smaller time 
separation having about half the statistical error. We therefore use tf/a = 8 or tf = 0.91 fm in what 
follows. 



4. Results 

We use the local electromagnetic current, j^ M = —^sy^s, which requires a renormalization 
factor Zy to be included. This renormalization constant is determined by the requirement that 
Geo at zero momentum transfer is equal to the charge of Q. in units of electric charge, that is 
Geo(0) = — 1. Our calculation, at this quark mass, yields a value of Zy = 0.727(1), which is 
reasonably close to the value obtained in Refs. ^] in the chiral limit. 

4.1 Electric charge form factor 

Our results for the electric charge form factor, Geo, are depicted in Fig. |5|(a). As can be seen, 
the momentum dependence of this form factor is described well by a dipole form 

Geo(Q 2 ) = - - * (4-1) 

with A| a fit parameter. In the non-relativistic limit the slope of the above dipole form evaluated 
at momentum transfer Q 2 = 0, is related to the electric charge root mean square (rms) radius 



(rlo) = -^G EQ (Q 2 ) 



(4.2) 

e 2 =o 
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Figure 2: The ratio 7?^ T (r,^,f) extracted for temporal source-sink separations tf/a = % and = 
10, using 50 gauge configurations. The results for tf/a = 10 are shifted to the left by one unit. 
We show results for current directions jx = 1,2,3 and low-momenta q: (0, 1,0) ^ and (1,0,0)^, 
respectively. 



From our dipole fit to the lattice data we determine A^o and obtain a value of (r| ) = —0.354(9) fm 2 



This is slightly greater in magnitude than the one reported in Ref. [jlOy, which was obtained in a 
quenched lattice QCD calculation (see Table [I]). Our value is expected to be higher since in a dy- 
namical lattice calculation meson-cloud effects are taken into account in addition to the fact that in 



Ref. [10] a heavy pion mass has been used. 



4.2 Magnetic dipole form factor 

Lattice results on the magnetic dipole Gm\ are shown in Fig. ||(b). As in the case of Geo, a 
dipole fit describes very well the Q 2 -dependence of Gm\- Fitting to the form 

G m {Q') = --^% (43) 

we can extract the anomalous magnetic moment of the Q. . By utilizing the lattice computed Q. 
mass, from Table [T], we calculate the magnetic moment in nuclear magnetons by using the relation 

Ma- = G M i (0) ^— = G m (0) — . (4.4) 

Our value of /i^- in nuclear magnetons is given in Table [l] and it is in agreement with the 
experimental value. It is also in agreement with two recent lattice calculations [|l(], [jl]]. The 
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Figure 3: (a) The electric charge form factor Geo at m n = 330 MeV. The line describes the dipole 



fit given in Eq. ( |4.1[ ). (b) The magnetic dipole form factor, Gmi, along with the dipole fit as given 
in Eq. (|3|) at m n = 330 MeV. The experimental datum for the magnetic dipole form factor, 
G m (0) = -3.60(8) GeV 2 [ffl], is also shown. 



calculation in Ref. JTo| ] is similar to ours in the sense that the three-point correlation function is 
calculated but the evaluation is carried out in the quenched theory and only at one value of Q 2 . 
In Ref. JTT| ] one employs a background field method to compute energy shifts using Nf = 2 + 1 
Clover fermions at pion mass of 366 MeV on an anisotropic lattice. 





Vol. 




[MeV] 


ma 

[GeV] 


Gm(0) 

[GeV 2 ] 




(4o) 

[fm 2 ] 


this work 
ref. [10] 
ref. |jll|] 


24 3 x 64 
20 3 x 40 
24 3 x 128 


200 
400 
213 


330 
697 
366 


1.763(21) 
1.732(12) 
1.650 


-3.58(10) 


-1.92(6) 
-1.697(65) 
-1.93(8) 


-0.354(9) 
-0.307(15) 


ref. \^ [PDG] 








1.672(45) 


-3.60(8) 


-2.02(5) 





Table 1: The Q. mass mn, the magnetic dipole form factor Gmi at Q 2 = 0, the magnetic moment Hq- and 
the electric rms charge radius (r| ). 



5. Summary 

Using appropriately constructed sequential sources the dominant Q. electromagnetic form 
factors Geo and Gmi are calculated with good accuracy even with a small sample of dynamical 
domain-wall fermion configurations. In the current calculation we neglected disconnected contti- 
butions. 
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The Q. magnetic moment is extracted by fitting the magnetic dipole form factor Gmi to a 
dipole form. We find a value that is in agreement with experiment |jl]]. The electric rms charge 
radius (r| ) is also computed and found to be larger than the value obtained in quenched QCD [|l0|]. 

We have also preliminary results on the subdominant electric quadrupole form factor Gei 
using the source of Eq. ( |3.10| ) but with our current statistics the errors are still large and no definite 
conclusion can be drawn. 

We will check for cut-off effects by performing the calculation of the form factors using dy- 
namical domain-wall fermion configurations at a finer lattice spacing. Although the light quark 
mass dependence is expected to be small, this needs to be checked at another, preferably smaller, 
value of the pion mass. 
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